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Abstract 

Two procedures for obtaining (extracting and constructing) the topological sig- 
nature of any multiply connected Robertson- Walker (RW) universe are presented. 
It is shown through computer-aided simulations that both approaches give rise to 
the same topological signature for a multiply connected flat RW universe. The 
strength of these approaches is illustrated by extracting the topological signatures 
of a flat (k = 0), an elliptic (k = 1), and a hyperbolic (/c = — 1) multiply connected 
RW universes. We also show how separated contributions of the covering isometries 
add up to form the topological signature of a RW flat universe. There emerges 
from our theoretical results and simulations that the topological signature arises (in 
the mean) even when there are just a few images for each object. It is also shown 
that the mean pair separation histogram technique works, and that it is a suitable 
approach for studying the topological signatures of RW universes as well as the role 
of non-translational isometries. 
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1 Introduction 



The astrophysical observations indicate that to a high degree of accuracy our universe is 
locally homogeneous and isotropic. Thus in the framework of the general relativity theory 
it can be described through a Robertson- Walker (RW) metric 

ds"^ = dt^ - R\t) da\ (1.1) 

where t is a cosmic time, and da'^ = dx^ + P{x) [ dO"^ + sin^ 9 dcj)^ ] with /(x) = X i sin x , 
sinh X , depending on the sign of the constant spatial curvature (/c = 0, ±1). Clearly R{t) 
is the scale factor that carries the unit of length. 

However, whether we live in a simply or multiply connected space and what is the 
size and the shape of the universe are open problems, which have received increasing 



attention over the past few years (see, for example, [|l| - and references therein). The 
most immediate consequence of multiply-connectedness of the universe is the possibility 
of observing multiple images of cosmic objects, whose existence can be perceived by the 
simple reasoning presented below. 

It is often assumed that the t = const spatial sections M of a RW spacetime manifold 
are one of the following simply connected spaces: Euclidean {k = 0), elliptic {k = 1), 
or the hyperbolic {k = —1), depending on the sign of the constant spatial curvature 
k. However, one can lift the restriction of simply-connectedness for the spatial sections 
by allowing M to be any one of the possible quotient manifolds M = M/T, where M 
stands for E^, or (depending on the sign of k), and F is a discrete subgroup of the 
full group of isometrics of the covering manifold M acting freely on M [^. The action 
of r tesselates M into identical cells or domains which are copies of what is known as 
fundamental polyhedron (FP). In forming the quotient manifolds M the essential point 
is that they are obtained from M by identifying points which are equivalent under the 
action of the discrete group F. Hence, each point on the quotient manifold M represents 
all the equivalent points on the covering manifold M. On the other hand, in the general 
relativity approach to cosmological modelling, the physicists assume that our universe can 
be modelled by a manifold, and thus a given cosmic object is described by a point p G M, 
which represents, when M is multiply connected, a set of equivalent points (images of p) 
on the covering manifold M. So, to figure out that multiple images of an object can indeed 
be observed if the universe is multiply-connected, consider that the observed universe is 
a ball Br^ C M whose radius Rh is the particle horizon, and denote by L the largest 
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lenght of the fundamental polyhedron P of M (P C M). Thus, when Rh > L/2, for 
example, the set of (multiple) images of an given object that lie in Br^ can in principle 
be observed. Obviously the observable images of an object constitute a finite subset of 
the set of all equivalent images of the object .Q 

General relativity (GR) relates the matter content of the universe to its geometry, and 
reciprocally the geometry constrains the dynamics of the matter content. As GR is a 
purely metrical (local) theory it cannot be used (without additional topological assump- 
tions as, e.g., asymptotical flatness or asymptotical locally anti-de Sitter |^) to settle 
the global structure (topology) of spacetime. One way to tackle the problems regarding 
the topology of the universe is through a suitable statistical analysis applied to catalogues 
of discrete cosmic sources to find out whether there are multiple images of cosmic objects, 
and eventually determine the shape and the size of the universe from the pattern of the 
repeated images. 

Cosmic crystallography (CC) is one of such statistical methods, which looks for 
distance correlations between the images of cosmic objects using pair separation his- 
tograms (PSH), i.e. graphs of the number of pairs of sources versus the distance between 
them. These correlations arise from the discrete isometrics of the covering group F, which 
give rise to the (observed) multiple images. The initial expectation behind the crystallo- 
graphic method was that these distance correlations would manifest as sharp peaks (called 
spikes) in PSH's, and that the spike spectrum (their positions and relative amplitudes) 
would be a definite signature of the topology. The first simulations performed for some 
specific flat manifolds appeared to confirm these expectations in that the PSH's 
corresponding to those particular manifolds presented distinct spike spectra. Afterwards, 



however, histograms were generated for the specific cases of Weeks [29| and one of the 



Best manifolds and no spikes were found. Their graphs show that, within the degree 
of accuracy of the corresponding plots, the PSH's of these hyperbolic manifolds exhibit 
no spikes. However, the spikes in PSH's can either be of topological origin or arise from 
purely statistical fluctuations. So, it is important to consider that the statistical noise 
in PSH's can both give rise to sharp peaks of purely statistical origin, and even hide the 
topological spikes, which are the sharp peaks that really matter in cosmic crystallography. 
Thus, the ultimate proof (or disproof) for the existence of the topological spikes in a class 
-"^For the precise and general condition for the existence of multiple images see Also note that for 
Rh > 2 rinj (here rinj is the injectivity radius, see e.g. j2l) for notation) multiple images can still arise. 
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of manifolds cannot rely only on graphs, since they do not constitute a proof but a clue 
or simply an indication.0 As a matter of fact, the proof for the existence of topological 
spikes ought to arise from a theoretical statistical analysis of the distance correlations to 
reveal the role in PSH's played by all types of isometrics of the covering group F. 

In considering discrete astrophysical sources in the context of multiply connected RW 
spacetimes, the observable universe is the region or part of the universal covering manifold 
M causally connected to an image of a given observer since the moment of matter- 
radiation decoupling. Clearly in the observable universe one has the set of observable 
images of the cosmic objects, denoted by O . A catalogue is a particular subset C C C, of 
observed images, since by several observational limitations one can hardly record all the 
images present in the observable universe. Our observational limitations can be formulated 
through selection rules which dictate how the subset C arises from O. Catalogues whose 
images obey the same well-behaved distribution and that follow the same selection rules 
are said to be comparable catalogues (for more details about the formalization of these 
concepts see p6|). 



In a recent article, in the context of Robertson- Walker spacetimes, Gomero et al. [26 



have derived an expression for the expected pair separation histogram (EPSH) for an 
ensemble of comparable catalogues C with approximately the same number of sources 
and corresponding to a manifold M. From the expression for the EPSH they have shown 
that the spikes of topological nature in PSH's are due to Clifford translations]^ whereas 
the other non-translational isometrics manifest as slight deformations of the EPSH corre- 
sponding to the underlying simply connected manifold. This general result holds regard- 
less of the 3-geometry of the spacelike sections of the RW spacetimes, and its restriction 
to the specific case of Euclidean and hyperbolic 3-geometries gives rise to two basic conse- 
quences: (i) that Euclidean manifolds which have the same translations in their covering 
groups exhibit the same spike spectrum of topological nature; (ii) that individual pair 
separation histograms (PSH) of hyperbolic 3-manifolds exhibit no spike of topological 
origin, since there are no Clifford translations in the hyperbolic geometry. From the ex- 
pression for the EPSH they have found it is also apparent that the contribution due to 
the multiply-connectedness (topological signature) must arise in PSH's even when there 

^For preciseness, it should be added that in two causes are presented for the absence of spikes in 
the simulated PSH corresponding to the Weeks manifold. Note, however, that the explanation in for 
the lack of spikes in PSH's does not seem to match with those presented in ref. |Q. 

Clifford translation is an isometry g of M such that for all p e M, the distance d{p, gp) is constant. 
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are only a few images for each object.^ These results have been formally derived from 
very general first principles in |2^, nevertheless they have not been thoroughly clarified 
through concrete simulations. 

In the quest for a means of reducing the statistical fluctuations well enough to clear 
up the signal of non-translational isometrics in PSH's, Gomero et al. |]2B[ also suggested 
the mean pair separation histogram (MPSH) as a very first approach to refine upon the 
crystallographic method. But, although the MPSH procedure rests on a well known result 



from elementary statistics, apart from the cases recently discussed in [^, no other 
MPSH has been explicitly built to show that the MPSH scheme is useful in simulations, 
and that it is a suitable approach for studying the role of non-translational isometrics as 
well as the topological signature of any 3-manifold of constant curvature. 

The most patent evidence of multiply-connectedness in PSH's is the presence of topo- 
logical spikes, which arise when the isometry is a Clifford translation. The other isome- 
trics, however, manifest as small deformations of the EPSH corresponding to the under- 
lying simply connected manifold. In computer-aided simulations, however, histograms 
contain statistical noise, which on the one hand can give rise to sharp peaks of purely 



statistical origin |3^, , on the other hand it can mask or even hide the tiny deformations 
due to non-translational isometrics. The most immediate approach to cope with these 
problems clearly is through the reduction of the noise in PSH's by using, for example, the 
MPSH scheme. Another possible way of facing them is by focusing in a more appropriate 
quantity rather than PSH's. Obviously one can also combine these approaches to tackle 
those problems. 



In this article, we firstly recast the theoretical results obtained in p6[ to show how 
one can extract the topological signature of any multiply connected manifold of constant 
curvature by using a suitable quantity, which turns out to be a constant factor times 
the difference ^exp{.Si) — ^^^^(sj) of the EPSH corresponding to the multiply connected 
manifold minus the EPSH of the underlying simply connected covering manifold. Secondly 
we show, through concrete simulations and based upon our theoretical results, two ways of 
obtaining (extracting and constructing) the topological signatures of multiply connected 
manifolds of constant curvature. There also emerges from our theoretical results and 
simulations how separated contributions of the covering isometrics are composed to give 



■^Note however, that whether or not the signatures can be distinguished from the statistical noise when 
studying a single observational catalogue is an independent question. 
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rise to the (complete) topological signature of a specific manifold. Further, we make clear 
that the topological signature arises (in the MPSH, i.e. in the mean) even when there are 
just a few images for each object. Finally, it is also illustrated that the MPSH procedure 
works by building MPSH's from simulated catalogues, and that it is a suitable approach 
for studying the topological signature as well as the role of non-translational isometrics. 

The scope of this paper is as follows. In the next section we set the notation, briefly 
recast the major result of ref. [^, and presents the theoretical grounds of the two ap- 
proaches for obtaining (extracting and constructing) the topological signatures of any 
multiply connected 3-manifolds of constant curvature. In section |^ we present and dis- 
cuss simulations that make apparent the following issues: (i) that the two procedures for 
obtaining the topological signature of multiply connected 3-manifolds discussed in sec- 
tion work, and give rise to the same signature; (ii) how the separated contributions of 
the covering isometrics add up to form the plain topological signature of RW universes; 
and (iii) that the topological signature arises in simulations (in the mean, i.e. through 
MPSH's) even when there are just a few images for each object. It is also shown in that 
section that the mean pair separation histogram technique works and that it is a suitable 
approach for studying the topological signature as well as the role of non-translational 
covering isometry of any multiply connected 3-manifold of constant curvature (spacelike 
sections of all RW universes). In section ^ we summarize our main results and briefly 
indicate possible approaches for further investigations. 



2 Topological signature 



In this section we shall first set the notation and briefly recast the major result of ref. ^6 
so as to show how one can extract a topological signature of any multiply connected 
manifold of constant curvature. 

Let us start by recalling that a catalogue C is a set of observed images, subset of the 
observable images, which are clearly contained the observable universe, which in turn is 
the part of the universal covering manifold M causally connected to an image of a given 
observer. 

Consider a catalogue C with n cosmic images and denote by //(s) the number of pairs 
of images whose separation is s. Consider also that our observed universe is a ball of 
radius a and divide the interval (0,2a] in m equal subintervals Jj of length Ss = 2a/ m. 
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Each of such subintervals has the form 



and is centred at 



5s 5s 

-^i = - Y > Si + yj ; « = 1,2, ...,m, (2.1) 

Si= {i- ^) 5s . 

The PSH is a normahzed function which counts the number of pair of images separated 
by a distance that hes in the subinterval Jj. Thus the function PSH is given by 

Hs,) = 1 E ^{s) ^J:v{s), (2.2) 

n{n i) OS ^gj. ^gj. 

and is clearly subjected to the normalizing condition 

m 

J2^isi)5s = l. (2.3) 

i=l 

In a multiply connected universe the periodic distribution of images on M due to 
the covering group gives rise to correlations in their positions, and these correlations 
can be couched in terms of correlations in distances between the pairs of images. The 
examination of the behaviour of these distance correlations can be made as follows. 

If one considers an ensemble of comparable catalogues with the same number n of 
images, and corresponding to the same 3-manifold M of constant curvature, one can 
compute probabihties and expected values of quantities which depend on the images in 
the catalogues of the ensemble. Clearly each catalogue of the ensemble gives a set of 
images distributed in the observed universe. In particular, one can compute the most 
important quantity for our purpose here, which is the expected number, rjexpisi), of pairs 
of cosmic images in a catalogue C of the ensemble with separations in Jj. This quantity is 
quite relevant because from it one has the expected (normalized) pair separation histogram 
which clearly is given by 

2 1 11 
n[n — 1) OS i\ 5s 

where obviously = n{n — l)/2 is the total number of pairs of cosmic images in C, and 
so the whole coefficient of rjexpisi) is nothing but a normalizing factor. 

If one denotes by Fg{si) and Fu{si), respectively, the probability that the images of a 

5f-pair[| and an uncorrelated pair be separated by a distance that lies in Jj, the expected 
^In line with the notation of | p6[ , when referring coUectively to correlated pairs we use the terminology 
F-pairs, reserving the name g-pair for a correlated pair corresponding to a specific isometry g Cz T. In 
other words, a g-pair is a pair of the form {p, gp) for any (fixed) isometry g. It should also be noted that 
the F-pairs are the same as the type II pairs of Uzan et al. . 



number rj^xpisi) can be expressed as 

Ve.p{si) = iV„ F^{si) + lT^Ng Fg{si) , (2.5) 

where F denotes the covering group F without the identity map, and where and Ng 
denote, respectively, the (total) expected number of uncorrected pairs and the (total) 
expected number of (7-pairs in a typical catalogue C of the ensemble.0 
Inserting eq. (|2.5|) into eq. (|2.4| ) one obtains 



^e.p{s.) = F^{s.) + i E f ^ • (2-6) 

9Gr 



However, from equation ( p.4| ) one is led to define the EPSH's corresponding to uncorre- 
lated pairs and associated to an isometry g, respectively, as 

^Ipis,) = i-lr/,Vs,) = lF„(sO, (2.7) 

= ^^r/e%(s,) = ^F,(s,), (2.8) 



where clearly = E., Vexpi^i) and Ng = Y.s, vLpi^i) ■ 

Thus using (^.7]) and (|2.8|), equation (^.61) reduces to either of the following forms: 



IN 

^exp{Si) = fin $Lp(Si) + E with /i« = ^ , = 2 ' ^^"^^ 

or equivalently 

<^>e.p(3i) = ^ [^.$:.p(3i) + E ^a^xpiSi)] With = 2^, z/, = ^ . (2.10) 

For simply connected manifolds M , since all pairs are uncorrelated, equation (|2.4| ) 
reduces to 

= ^ ^ = ^ ^.c(s.) , (2.11) 

where Fsdsi) is the probability that the two objects in M be separated by a distance that 
lies in Jj.Q 

Now since the pairs of cosmic images are either correlated (F-pairs) or uncorrelated 
we must have 

Nu + l^Ng = N , (2.12) 
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For a formal proof that the decomposition ( ) can always be made see 



^It should be noticed that although all pairs in the catalogues corresponding to simply connected 
manifolds are uncorrelated the probabilities Fu{si) and Fsc{si) are not the same. 
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which on the one hand leads to 



liu + T,lig = l, (2.13) 

and on the other hand it gives rise to an alternative expression for the EPSH ^exp{^i) iii 
terms of ^l%p{si) . Indeed, using (p^) and ( p.l3|) one easily obtains 



^exp{Si) = ^f^p{Si) + /i, [ Kxp{Si) - Klp{Si) ] + E yt^g [ " Klp{Si) ] ■ (2.14) 

From equation ( ^.141 ) one can obtain an expression for the topological signature 
ip^{si) = (n — l) l^expisi) — ^llp{si)] . Indeed, using ( p^.l4D together with the definition of 
fiu and one easily obtains 



where 



and 



^'(s,) = (n - 1) [ $,.p(s,) - $::,(s,) ] = ¥.^(s,) + /(s,) , (2.15) 

^''{s,) = i^j^:,pis,)-^::pis,)] , (2.16) 



^^is,)= Y.^Anxpis.)-K:pis.)], (2.17) 
and where we have used the definitions of and Ug given by ( p.lO| ). 



An important word of clarification is in order here: although throughout this article 
we loosely use the terminology topological signature of a manifold and of a universe, it 
should be stressed that the topological signature f^{si) depends on the pair 3-manifold 
and observed universe Ba , as well as on the relative position of the FP with respect to 

To show how one can extract through computer-aided simulations the topological 
signature (p^{si) = {n — 1) [^exp{si) — ^llp{si)] of multiply connected universes, a quite 
important point to be noted is that the EPSH is essentially a typical PSH from which the 
statistical noise has been withdrawn. Hence we have 

<^>{Si) = <l>exp{Si) + p{Si) , (2.18) 

where $(sj) is a typical PSH constructed from C and p{si) represents the statistical 
fluctuation that arises in the PSH $(si). Alternatively from ( ^IsD one has that a PSH 
is an EPSH plus the statistical fluctuation. 

9 



Using now the decomposition ( 2.18 ) together with ( p.l5 ) one easily obtains 



{n - 1) [ $(5,) - <^-'%s,) ] = if^Si) - (n - 1) [p{si) - p^^(s,) ] , (2.19) 

where p(sj) and p^^{si) are the statistical noises which arise in the PSH's corresponding 
to the multiply connected and the covering manifolds, respectively. 

Clearly the right-hand side of ( ^.19| ) gives the topological signature ^p'^{si) intermixed 
with the two statistical fluctuations p(sj) and p^'^{si). In practice, however, one can ob- 
viously approach to the suitable quantity [^exp{si) — $g^p(sj)] by reducing the statistical 
fluctuations, through any suitable statistical method to lower the noises. The simplest way 
to accomplish this is to use several comparable catalogues, with approximately equal num- 
ber of cosmic images, for the construction of a mean pair separation histogram (MPSH). 
For suppose we have K computer-generated catalogues Ck {k = 1,2, K) whose PSH's 
for a given value of m are given by 

ntiuk - 1) 5s 

where is the number of images in the catalogue Ck- The MPSH defined by 

<^s,)> = ^j2Ms^) (2.21) 
^ k=l 

is such that in the limit K ^ oo one has ^exp{si) — < ^{si) > ■ Elementary statistics 
tells us that the statistical fluctuations in the MPSH are reduced by a factor of 1/^/K, 
which makes clear that MPSH is a suitable approach to deal with the above-mentioned 
fluctuations. 

In brief, the use of the MPSH to extract the topological signature f^{si) consists in 
the use of K (say) computer-generated comparable catalogues (with approximately the 
same number n of images and corresponding to the same manifold M) to obtain the mean 
pair separation histogram < $(sj) > and analogously to have < $*'^(sj) > ; and use them 
as an approximation for ^expi^i) and $e^.p(sj) , to construct the topological signature 
ip^{si) ^ {n — 1) [< $(si) > — < $'*^(sj) >] . Clearly the greater is the number K of 
catalogues the lower are the statistical noises p(sj) , p*'^(sj) and obviously the better are 
the approximations < $(sj) > ~ ^expisi) and < ^^'^{si) > ^ ^Hpisi) ■ 

An improvement of the above procedure to extract the topological signature ^p^{si) 
comes out for the cases one can derive the expression for the PSH's $g^p(sj) corresponding 
to the simply connected covering manifolds. In a recent work (see also ||3^) the explicit 
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formulae for the expected pair separation histogram functions ^llp{si) corresponding to 
an uniform distribution of objects in the covering manifolds E^, and 5*^ have been 
obtained.0 Thus, for such universes with homogeneous distribution of objects one has from 
the very beginning p'^^(sj) = , and obviously the topological signature can be rewritten 
as if'^isi) ~ - 1) [ < > - ^Zj,is,) ] . 

For the sake of selftcontainedness we present below for use in section ^ the explicit 
expressions for the EPSH's found in refs. [0, ^ corresponding to simply connected 



Euclidean, hyperbolic and elliptic universes with a uniform distribution of cosmic objects. 
Euclidean Universes 



^llpi^^ = ll^ - sf {s + 4a) , (2.22) 



which holds for s G (0, 2a] . 
Hyperbolic Universes 



S sin fi ^ s 

= (sinh2a-2a)^ sech(./2) sinh(a - s/2) - (a - s/2) ] , (2.23) 



which also holds for s G (0, 2a]. 
Elliptic Universes 



i sm^ s 



*eL(«'^) = 77^ ^— -r {2a-sin2a-7r + e(27r-2a-s) [sin2a + TT 

^""P (2a - sm2a)2 ' ^ 

-a -s/2 -cos a sec(s/2) sin(a - s/2) ] } , (2.24) 

which holds for all a G (0, tt] and s G (0, min(2a, vr)] , and where denotes the well-known 
Heaviside function. 

It should be noticed that in ( ^.24| ) as well as in ( f^.22D and ( [^.23D we have made explicit 



that the EPSH's depend upon the radius a of the observed universe Ba ■ 

From equations ( |2.15| ) - ( p.lTj ) it is also clear that another approach to obtain (con- 
structive approach) the topological signature of a multiply connected manifold is by 
considering the sum of the terms on the right-hand side of (|2.15|) , namely V'^(si) = 

J^ui^expisi) - ^Zpisi)] and ip^{s,) = Eg^fJygl^Lpisi) - ^Zp{s^)]■ In the next section 
^Actually, they have found the explicit expressions for the EPSH corresponding to each (Euchdean, 
hyperbohc and eUiptic) observed (homogeneous) universe, which is a spherical ball Ba with radius a 
fulfilled with an uniform distribution of cosmic objects. 
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we shall perform simulations to obtain (construct) the topological signature of a given 
manifold using also this alternative procedure and make clear that one ends up with the 
same signature obtainable through the first approach, as it should be expected from the 
outset. In doing so we will be showing in addition that indeed (|2.15| ), with (p.l6|) and 
( |2.17D , holds within the accuracy of our plot, of course. 

To close this section it should be noted that in a recent work, in studying PSH's cor- 
responding to one of the Best manifold, Fagundes and Gausmann have considered 
the similar quantity — <l>'*^(si), instead of the quantity ^exp{si) — we are 

interested here. However, apart from the fact that they have not presented any theo- 
retical argument sustaining their suggestion, from (|2.19|) one has that their scheme gives 
rise to a fraction l/{n — 1) of the topological signature f^{si) plus (algebraically) the 
fluctuations corresponding to both PSH's involved, namely that which arises from the 
multiply connected Best manifold and that originated from the covering space. 



3 Topological signature from computer-aided simu- 
lations 

In this section we shall present and analyze simulations that clarify the following 
points: (i) how one can obtain (extract and/or construct) the topological signature of 
multiply connected manifolds; (ii) how separated terms which arise from the covering 
isometries are composed to form the topological signature of a manifold; and (iii) show 
that the topological signature arises in simulations (in the mean) even when there are 
just a few images for each object. In the clarification of theses points we shall use the 
MPSH technique to make apparent that it is a suitable approach to study the topological 



signature and the role of non-translational isometries as it has been stated in PBI and 



illustrated in ||3T| for the case of flat manifolds. 

The first series of computer-aided simulations concerns a compact orientable Euclidean 
manifold of class Qq in Wolf's classification [^. Recently it has been found a set of 
expressions for the face-pair generators of a fundamental polyhedron (FP) of a cubic 
manifold of this class. We shall denote this Euclidean manifold by 7^ in agreement with the 



notation used in ||36| and , wherein the cubic FP and the pairwise faces identification 
are shown in figure 1. Relative to a coordinate system whose origin coincides with the 
center of the FP, the actions of the generators a, (3 and 5 on a generic point p = {x, y, z) 
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were shown to be described by 



25. 



ap = {x + L, —y, —z) 
PP = {-x, z + L, y) , 
Sp = z,y + L), 



(3.1) 
(3.2) 
(3.3) 



where L is the edge of the cubic FP. Clearly the actions of the inverses of these generators 
are given by 



In the simulations corresponding to the cubic manifold the centre of the FP was 
taken to be the origin of the coordinate system, and to coincide with the centre of the 
observed universe Ba , whose diameter 2a is such that 2a = L \f2 ~ 1.41 L . It should be 
noted that with this ratio for ajL and for s G (0,2a) one has only the contribution of 
non-translational isometrics for the topological signature. Indeed, according to eqs. ( p.l|) - 
( |3.6|) the translations of shortest length are due to and their inverses, and are 

translations of L \/2 units. So, as we shall take < s < 2a no spike will appear in our 
PSH's for this cubic manifold.0 

To show how one can extract the topological signature [n — 1) [$ezp(si) — '^'exp("^«)] — 
<f^{si) of 74 , as well as to illustrate the MPSH procedure, and to make clear that the 
topological signature arises in simulations when there are a few images of some cosmic 
objects, we have written a program whose input are the number K of catalogues, the radius 
a of the observed universe Ba , the number m of subinterval (bins), and the number Ug 
of objects inside the FP (seeds). The program generates K different catalogues, starting 
(each) from the same number n<j of homogeneously distributed seeds inside the FP. For 
each bin Jj of width 6s = 2a/m it counts the normalized number of pairs I]%(s) for 
all catalogues k from 1 to K. Finally, it calculates the normalized average numbers of 
^Note, en passant, that the topological signature depends on the pair: 3-manifold and observed 
universe Ba , and also on the relative position of the FP with respect to observed universe Ba ■ So, had 
we considered, for example, a ball with double diameter Aa — 2 V2 ~ 2.83, and kept the same relative 
position of FP and Ba , the topological signature ip^{si) would exhibit spikes at s = \/2 ~ 1.41, s = 2, 
and s — V6 ~ 2.45 and so forth (see ref. |3^| for more details). 



a p 



{x - L, -y, -z) 



(3.4) 



(-x, z,y -L) , 
(-X, z-L,y). 



(3.5) 
(3.6) 
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pairs for all Sj G (0, 2a) , finding therefore, according to ( 2.21 ), the mean pair separations 
histogram <$(si)> over i^' catalogues. 

Now, on the one hand we know that the mean pair separation histogram < $(sj) > is 
approximately equal to the expected pair separation histogram ^exp{si) for large enough 
number K. On the other hand, we have that the explicit expression for ^llp{si) corre- 
sponding to a uniform distribution of objects in a Euclidean universe is given by eq. ( |2.22| ). 
Thus, the plot of (n — 1) [ < > — $g^p(sj) ] gives, up to a statistical fluctuation p{si), 
a topological signature of the 74. As a matter of fact, since not all catalogues generated 
from a fixed number of seeds inside the FP have the same number n of images, in our 
program instead of n we have used the average number <n> of images per catalogue to 
plot the topological signature v^'^(si) . 

Using the above-described program we have performed simulations for the manifold 
74 with L = 1 in an observed universe with radius a = V2/2 ^ 0.71 , 6s = 0.01, and with 
different number of seed objects uniformly distributed in the FP. Figure la is the graph 
of the topological signature {<n> — 1 ) [< > — ^^^^(sj)] for the manifold 7^ and 

was obtained using the MPSH procedure for K = 16000 catalogues, and Ug = 15, which 
corresponds to an average number of images per catalogue <n> ^ 23. Figure lb shows a 
graph of the topological signature for the same manifold and obtained through the MPSH 
scheme for identical number of catalogues, but now the number of seeds was taken to be 
Hg = 100, which corresponds to < n > ~ 153 images. These figures clarify the following 
relevant points: First, that the MPSH approach is indeed a suitable approach to reveal the 
topological signatures of non-translational isometries; Second, as long as one can extract 
the topological signature for a small numbers of seeds (n^ = 15) and it is essentially the 
same obtained for a fairly large number of seeds {ug = 100), it becomes clear from these 
concrete simulations that the plain topological signature arises in simulations where there 
are just a few images for each object. 

According to our earlier discussion, from equation ( p.l5| ) - (|2.17| ) it also is clear that 
an alternative approach to obtain (constructive approach) the topological signature of 
a multiply-connected manifold is to consider the sum of the terms on the right-hand 
side of ([2J[^ ), namely ip^{si) = z/, [<l>^,p(s,) - ^llp{s,)] and ip^{s,) = Eg^fJ^gi^Lpisi) - 
^expi^i) ]• have also performed simulations following this approach to determine the 
topological signature of the manifold T4 again with L = 1 and a = ^/2/2 ~ 0.71 . In 
doing so we intend, on the one hand, to numerically verify equation ( p.l5| ), on the other 
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hand, to illustrate how the separated terms which arise from the covering isometrics sum 
up to give the topological signature. In what follows we shall report more details of the 
simulations involved in this second approach. 

To make explicit the contribution $f^p(sj) of a specific isometry g = a (say) to the 
topological signature term (p^{si) [see eqs. ( ^.8| ) and ( |2.15| ) - ( ^.171 )] in computer-generated 
catalogues we have written a program whose inputs are the total number A^^ of correlated 
pairs, the length L, the radius a of the universe Ba and the number m of subintervals 
(bins); and starts by setting a counting variable k = and has the following steps: 

(i) randomly and homogeneously take an object p (say) inside FP 

(ii) apply the isometry a and check if the image ap is still inside the ball Ba ■ If no, 
return to step (i), if yes let k = k + 1 and go to step 

(iii) calculate the distance s between p and ap, then 

(iv) call a procedure that computes for all s G (0, 2a) the normahzed sum of pairs v{s) 
whose images are separated by a distance that lies in a bin Jj with width 6s = 2a/ m; 

(v) check if the k is equal to Ng. If no, return to step (i); otherwise plot the sequence 
of pairs [x = {j — 1/2) 6s , y = J2v{si)/Ng] for j = 1, . . . ,m and stop. 

We have used this program to compute the contribution of the isometrics a and (3 for 
Ng = 16000, m = 142 {6s = 0.01), a = ^2/2 and L = 1. Figures 2a and 2b show, 
respectively, the graphs of < $°'(sj) > and < $^(sj) > for the isometrics a and /3, whose 
actions are given, respectively, by ( ^.1| ) and ( ^.2| ). The fact that the action of isometry a, 
for example, gives rise to g-pairs whose minimum separation is 1 can be easily understood 
from the expression for the distance d between p and a p, namely from d{p, ap) . Indeed 
it is trivial to show that d{p,ap) > 1. Similarly, form ( |3.2| ) one clearly has d{p, jSp) > 
\/2/2 ^ 0.71 , making clear that the /3-pairs have also a minimum separation.]^ 

For the ratio a/L = V2/2 and relative position of FP and Ba that we have chosen, only 
the isometrics given by ( p.l|) - (|3.6|) will contribute to the term (f^{si). Furthermore, from 

Actually, it can be explicitly shown that the set of starting points at which the non-translational 
isometrics begin to contribute for the topological signature are given by cf — 0.5, 1,4.5,9, 12.5, and so 
forth [^. Thus, it is clear that the inclusion or exclusion of the contribution of a specific isometry 
depends on the radius a of the observed universe Ba , reinforcing the fact that the topological signature 
depends upon the ball Ba ■ 
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these equations one clearly has $exp (sj) = $g^p(sj) and (sj) = $f^p(si). Again, due to 
the symmetry of eqs. ( |3.1|) - (|3.6| ) one easily obtains that $f^p(si) = = $f^p(si) . 



Therefore, the second term on the right hand side of ( p. 151) reduces to 



/(s.) ^ i/,[2 <<l."(s.)> +4 <$^(.,)> -6$::,(s,)] . (3.7) 

Finally, for a/L and the relative position of FP and Ba we have chosen, according to 
equations ( p.l| ) - ( p.6| ) one easily has that for an uniform distribution of objects the 
coefficients Vg are equal for all isometrics g = a, (3,6 and their inverses. The computation 
of the value Ug reduces to the calculation of quotient of volumes, and turns out to be 
Ug ^ 0.1161165234 (for more details on how one can calculate Ug see [0). Figure 3 
shows the graph of $g^p(si) for a = a/2/2 as given by ( |2.22| ). Figure 4a shows the graphs 



of this (f^{si), given by (p.7|) , for 5 = 0.01, a = v^/2 and L = 1. It should be stressed 
that this figure is nothing but the above-described suitable combination of figures 2 and 3 
multiplied by the constant factor Ug . 

Again for the same relative position between FP and Ba , the contribution of the term 
^p^{si) to the topological signature [ip^{si) ] of 7^ in shown in figure 4b, and was obtained 
for the following set of inputs: m = 142, a = V2/2, L = 1, and K = 16000. 

Figure 5 gives the sum of the contributions due to ^p^{si) and v^^(sj) . The comparison 
between figures la and lb with figure 5 shows on the one hand that both procedures to 
obtain the topological signatures work; on the other hand these figures constitute to a 
certain extent a numerical check for the expression ( |2.15| ). Furthermore, figures 2 to 5 



also show how the separated terms, which arise from the covering isometrics, give rise 
to the topological signature looked upon as the sum of the terms in the right-hand side 
of (|2T5|) . 



To make explicit that the above procedure employed to extract the topological signa- 
ture of the Euclidean {k = 0) manifold T4 can be similarly applied to any other classes 
of 3-manifolds of constant curvature, we have performed computer-aided simulations for 
an elliptic {k = 1) as well as a hyperbolic {k = —1) manifolds. In the remainder of this 
section we shall briefly report the results of our simulations without going into details of 
the calculations (and programs) for the sake of brevity. 

We recall, firstly, that all (locally homogeneous) manifolds with k = 1 are known 
and that they are the simply connected compact covering manifold M = S^, the (compact) 
projective space = / Z2 and the following (compact) quotient manifolds M = S^/T, 
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where T is one of the following groups: (i) the cyclic group of Zp of order p > 2; (ii) 
the binary dihedral groups D* of order 4r (r > 2); (iii) the binary polyhedral groups: T* 
(where T is the symmetry group of the regular tetrahedron), O* (where O is the symmetry 
group of the regular octahedron), and J* (where / is the symmetry group of the regular 
icosahedron) ; (iv) the groups of the form D* x Zp with p > 2 and r > 2, although for 
certain values of r, D* x Zp can act on in two different ways; and (v) the groups of 
the form H x Zp {p > 2), where H = T*,0*,I*,T*, and the groups T* are non-cyclic 
subgroups of T*. 

We have performed computer simulations for the specific elliptic 3-manifold S'^/Z^, 
whose volume 2 7r^i?^/5 is one fifth of the volume of S'^ . A FP (tetrahedron) together 
with the pairwise faces identifications is given by Weeks [Q. We have taken as the 
observed universe the whole covering space S*^, i.e. a solid sphere with radius a = tt 
[R = 1 in ( |1.1|) ]. Thus all catalogues in our simulations for this manifold have the same 
number of images. Figure 6 shows the graph of the topological signature '{>^{si) = 
{n — 1)[< $(sj) > — ^llpisi)] for a manifold S'^/Z^ with the edge of the tetrahedron 
L ^ 1.82 and for m = 180, n = 100 images (n^ = 20), K = 3000 catalogues, and where 
the expression ( p.24| ) for ^llp{si) has been used. Clearly the existence of topological 
spikes in PSH gives rise to spikes in ip^{si). Now since the observed universe is the whole 
unitary sphere and the graph of ^p'^{si) shown in figure 6 presents no spike, it becomes 
apparent that the covering group F of this 3-manifold S^/Z^ has no translations. 

Contrarily to the case of elliptic space-forms a complete classification of the 3- 
dimensional hyperbolic manifolds has not yet been accomplished. However, from a topo- 
logical viewpoint the negatively curved (hyperbolic) universes are generic in the sense that 
most 3-dimensional manifolds can be viewed as homogeneous negatively curved and com- 
pact The volume of the covering manifold is infinite, but the quotient manifolds 
H^/r can either be finite of infinite. 

We have performed computer simulations for the specific compact hyperbolic 3-man- 
ifold known as Seifert- Weber dodecahedral space, which is obtained by identifying or 
glueing the opposite pentagonal faces of a dodecahedron after a rotation of 3 vr/S . Figure 7 
shows the graph of the topological signature (p^{si) = {<n> — 1) [< $(sj) > — ^l'^p{si)] 
for this hyperbolic space where the centre of the dodecahedron was taken to coincide with 
the centre of the observed universe Ba , whose diameter is 2a = 2.88. The length L of 
the edges of the pentagonal faces and the height H of the dodecahedron are such that 
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L = H = 1.992769 , where the lengths are measured with the hyperbohc geometry (|1.1| ) 
with R = 1. We have taken m = 100 bins, Ug = 10 seeds, (<r2>~ 18) K = 16000 
catalogues, and used the expression ( p.23|) for ^l'!^p{si). Since the existence of topological 
spikes in PSH gives rise to spikes in ip^{si) and reciprocally, figure 7 also shows that the 
PSH of even this highly symmetrical hyperbolic manifold presents no topological spikes, 
clarifying the conjecture made in ||3^ and in agreement with the results of ||26[ . 

It should be stressed that the unit of lengths used in all simulations of this section 
could certainly have been taken to be of hundreds or thousands of Mpc obviously without 
changing the patterns of the simulations we have performed. We have avoided large 
numbers for the sake of simplicity. 

To close this section we mention, for the sake of completeness, that small multiply- 
connected flat models (like those with 7^ topology) may be consistent with COBE CMB 
constraints (this is still a controversial point, see in this regards, e.g.. Levin et al. §], 
Roukema and Inoue |^). Nevertheless, according to the current values of the cos- 
mological parameters, a spherical universe is unlikely to be completely observable. Indeed, 
from the redshift-distance relation for Friedmann-Lemaitre-Robertson- Walker with dust 
and cosmological constant, one can easily obtain that a spherical universe with matter 
current density parameter VLmo = 0.3 would require that Q^o > 1-35 for the antipodal 
point to be visible at a redshift below z = 1000, i.e. for 5*^ to be inside the sphere 
representing the microwave background. Similarly with Q^o = 0-3 would require that 
Q\Q > 1.62 for an antipodal point to be visible at a redshift below z = 3, i.e. for to 
be coverable by galaxies and quasars. These range of values for Q/^q are in conflict with 
current constraints on Q^o and on = Qmo + ^Ao (see, for example, Lange et al. [|^, 



Bond et al. 44 , Balbi et al. [45 1 and Bernardis et al. [Eq 



4 Concluding Remarks 

If we live in a multiply connected Friedmann-Lemaitre-Robertson- Walker (FLRW) uni- 
verse the sky may show multiple correlated images of cosmic objects. The image correla- 
tions are dictated by the discrete isometrics of the covering group P of the 3-manifold used 
to model its space sections. The periodic distribution of images gives rise to correlations 



in their positions. In the crystallographic method [2S] these correlations are couched in 



18 



terms of correlations in distances between the images. Actually, the method of cosmic 
crystallography (CC) looks for distance correlations between cosmic images using pair 
separations histograms (PSH), with normalized function given by ( p. 21) , and whose graph 
gives the (normalized) number of pairs of sources versus the distance between them. 

The primary expectations in CC were that the distance correlations would manifest 
as topological spikes in PSH's, and that the spike spectrum would be a definite signature 
of the topology. Although the first simulations performed for specific flat manifolds 



appeared to confirm the initial expectations [2S], histograms generated subsequently for 



hyperbolic manifolds p9| , |30|| revealed that the PSH's of those manifolds exhibit no spikes. 
Concomitantly, a theoretical statistical analysis of the distance correlations in the PSH's 
was performed, and a formal proof was presented that the spikes of topological origin in 
PSH's are due to translations alone This result explains the absence of spikes in 

the PSH's of hyperbolic manifolds, and also gives rise to the fact that Euclidean distinct 
manifolds which admit the same translations on their covering group present the same 
topological spike spectrum |3l| . 

Although the set of topological spikes in PSH's is not definite topological signature 



and is not sufficient for distinguishing even between some compact flat manifolds [pi 



the most striking evidence of multiply-connectedness in PSH's is indeed the presence of 
topological spikes, which arise only when the isometry is a Clifford translation. The 
other isometrics, however, manifest as tiny deformations of the expected pair separa- 
tion histogram (EPSH) corresponding to the underlying simply connected universe. In 
computer-aided simulations, however, histograms contain statistical fluctuations, which 
can give rise to sharp peaks of statistical origin, or can hide (or mask) the tiny deforma- 
tions due to non-translational isometrics. 

The most immediate approach to cope with fluctuation problems is through the reduc- 
tion of the noise in PSH's by using the mean pair separation histogram (MPSH) scheme to 
obtain <$(si) > rather than a single PSH $(si) . However, in most of the computer sim- 
ulations we have performed, for a reasonable number of images (ra ^ 60 , or so) and when 
there is no topological spikes (no translations), the graphs of the ^exp{si) ^ < $(sj) > 
and $g^,p(sj) are essentially the same (the deformation of the EPSH due to the non- 
translational isometrics are indeed rather tiny!), making clear that in practice ^exp{si) is 
not a suitable quantity for revealing the topology of multiply connected universes. 

In this work we have studied two ways of obtaining (extract and construct) the 
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topological signature of any multiply connected RW universe. The important points 
in the first approach are: (i) the use of (|2.14|) to introduce a new quantity ip^{si) = 
{n — 1) [$exp(si) — ^^expi^i)] ) "which turns out to be suitable for revealing the topologi- 
cal signature; (ii) the supplementary use the MPSH technique to drastically reduce the 
statistical noises and therefore improve the approximations < > ^ ^exp{si) and 

The second way of obtaining the topological signature of a multiply connected RW uni- 
verse (constructive approach) is based on the explicit expression for (p{si) given by ( |2.15| ) - 
( p.l7|) . In this approach the topological signature is obtained by considering the sum of 
the terms on the right-hand side of (|2.15|) together with the MPSH technique to reduce 
the statistical noises. 

We have also shown through concrete computer-aided simulations and based upon our 
theoretical results, that the two ways of obtaining the topological signatures of multiply 
connected RW manifolds give rise to the same topological signature, as one would have 
expected from the outset. Furthermore, the strength of these approaches has been shown 
by extracting the topological signatures of a flat (A; = 0), an elliptic (/c = 1), and a hy- 
perbolic {k = —1) multiply connected RW universes. There emerges from our theoretical 
results and simulations that the topological signature arises (in the mean) even when 
there are just a few images for each object. 

It should be emphasized that although in the simulations of section |^ we have assumed 
that the matter is homogeneously distributed in the universe, the technique we have dis- 
cussed in this article can be similarly used for other types of matter distribution. In the 
more realistic cases in which the galaxies (or cluster of galaxies) cluster, the matter distri- 
bution is clearly inhomogeneous. This inhomogeneity has been cast in terms of statistical 
indicators such as correlation functions, which can be used either for the computation 
of the expression for $g^p(sj) or to simulate the clustered distribution of matter in the 
simply-connected case to have the approximate contribution ~ <l>g^p(sj). Once either of 
these steps has been carried out one can proceed along the lines described in this paper. 

It should be stressed that the ultimate step in most of the statistical approaches 
to extract the topological signature is the comparison of the signature obtained from 
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As a matter of fact, although in general necessary, in the present article we did not use the approx- 
imation < $'*'=(si) >~ ^l'^p{si) , but the explicit expressions for ^l'^p{si) given by ( 2.22 ) - ( 2.24 ) for 
universes fulfilled with an uniform distribution of cosmic objects. 
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simulated catalogues against similar ones generated from real catalogues. To do so one 
clearly has to have the simulated patterns of the topological signatures of the manifolds, 
which can be achieved by the approaches discussed in this paper. Note, however, that the 
MPSH technique (needed in both approaches) is restricted to simulated catalogues, since 
it is impossible in practice to have a reasonable ensemble of comparable catalogues of 
real cosmic sources to calculate mean and expected values of the relevant quantities. An 
approach to face the important remaining fluctuations problem is to study quantitatively 
the noise which arises in ( |2.19|) in order to develop suitable filters for those statistical 



noises which naturally arise in PSH's built from real catalogues. 

Note, however, that we have assumed that all cosmic objects of our interest are point- 
like and have long lifetimes so that none was born or dead since the time corresponding 
to the redshift cutoff of the catalogue. Moreover, we have also assumed that all objects 
are comoving, so that their worldlines have constant spatial coordinates. Although sim- 
pliflying these assumptions are commonly used in the literature and are very useful to 
study the observational consequences of a non-trivial topology for the universe. 

The main sources of uncertainties of the known statistical approaches to determine 
the topology of our universe from discrete sources are (see, for example, and [^] ): 



(i) the cosmological parameters flmo and ^ao are not accurately known or determined, 
so one cannot accurately compute distances from redshifts. The uncertainty in 
the determination of the cosmological parameters give rises to uncertainty in the 
determinations of the positons of the sources, which in turn lead to errors in the 
determination of distances between pairs of sources. 

(ii) the cosmic objects are excactly not comoving, so multiple images are not where they 
ought to be: in real catalogues there are uncertainties in the positions of the objects 
due to their peculiar velocities. 

(iii) most cosmic objects do not have very long lifetimes, so there may not exist images of 
the same object at very large different distances from one of our images (observer). 

(iv) in most studies of such statistical approaches to cosmic topology the catalogues are 
taken to be complete. Real catalogues, however, are incomplete: objects are missing 
due to selection rules, and also most surveys are not full sky coverage surveys. 

It should be noticed that if the universe turns out to be small with a covering group 
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containing Clifford translations, the collecting correlated pairs method (CCP method), 
which has been recently devised by Uzan, Lehoucq and Luminet (see also |^), can 
be used in conjunction with the cosmic crystallographic method to determine both the 
cosmological parameters flmo and f^Ao as well as the topology of the universe. In the 
remainder of this paper we shall discuss this point in more details. 

Instead of grouping the pairs of sources according to their separations to extract the 
signature of the correlated pairs, in the CCP method one uses the basic feature of the 
isometrics, i.e., that they preserve the distances between pairs of images regardless of 
whether they are correlated or not. Thus, if {p, q) is a pair of arbitrary images (correlated 
or not) in a given catalogue C, then for each g & T such that the pair {gp^ gq) is also in C 
we obviously haveQ 

d{p,q) = d{gp,gq) . (4.1) 

This means that for a given (arbitrary) pair {p, q) of images in C, if there are n isometrics 
(7 G r such that both images gp and gq are still in C, then the distance d{p, q) will occur 
n times. 

The CCP method consists essentially in a count of the number of times each pair 
distance between images is repeated in a given catalogue.[^ Actually, the CCP-index TZ is 
defined to be the quotient between the result of this count and the total number of pair 
of sources in C minus 1 , i.e. (A^ — 1) in the notation of the second Section. 

Now, given that in a catalogue corresponding to a simply connected universe the 
probability that a given pair separation be repeated is zero, a non-zero CCP-index is an 
indicator of a non-trivial topology (multiply-connectedness) of our Universe. However, as 
it has been discussed in [^, the CCP-index TZ gives no hint on what is the topology of 
the spatial sections of the universe. 

It should be noticed that cosmic crystallography, its variant approach for simulated 
catalogues discussed in this work, as well as the CCP method, are them all based on the 
distance determination between cosmic images, and therefore they are rather sensitive 
to the (precise) values of the cosmological parameters f2mo and i^Ao. However, as it has 



been discussed in ref. |47|, in a small universe the graph of TZ in terms of the cosmological 



"'^^The pairs for which eq. (4.1) holds have been referred to by Uzan et al. |47| as type I pairs. 

^^Note that one necessarily has a nonnuU CCP-index when there are pattern repetitions in a catalogue 

C as that of the scheme used by Roukema [^] in the search for repetition of local quintuplets of quasars in 

our universe. As he has pointed out his choice for quintuplet rather than any n-tuples was a compromise 

between the data and computing power available, but the isometrics clearly preserve the local n-tuplets. 
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parameters and Oao should exhibit a resonance spike at the correct values of these 
parameters. Therefore, one can use the CCP method to achieve a precise determination 
of these parameters. So, if the Universe turns out to be (small) Euchdean or spherical, for 
example, one can use the values of Q,rno and JIao determined through the CCP method, 
and then apply cosmic crystallography (CC) to determine the spike spectrum in the PSH. 
If the spike spectrum does not uniquely determine the topology of our Universe (as in 
the case for compact Euclidean manifolds), one can use again the CCP-index in order to 
distinguish among the possible candidates. Indeed, since the CCP-index is essentially a 
multiple count of the number of isometrics that give rise to repeated distances of images 
in the observed universe, then the CCP-index corresponding to a given covering group is 
larger than the one corresponding to any of its subgroups. Therefore, detailed simulations 
of the CCP method apphed to each hkely candidate (manifold) can eventually permit to 
single out one of the possible manifolds (candidates) , determining therefore the topology 
of the universe. This combined method (CCP plus CC and/or its variant approach for 
simulated catalogues) is currently under investigation, and our results will be shortly 
published elsewhere. 

Captions for the figures 

Figure 1. The topological signature (p^{si) = (< n > — l)[<<l>(si) > — ^Ixpi^i)] (ob- 
tained through the MPSH technique) of the cubic manifold with L — 1, in an 
observed universe with radius a ~ 0.71 . The horizontal axis gives the pair separa- 
tion s while the vertical axis furnishes the normalized number of pairs. In (a) the 
number of seeds is = 15 and corresponds to an average number of images per 
catalogue < n > ~ 23 . In (b) the number of seeds is Ug = 100 and corresponds to 
an average number of images per catalogue <n> ~ 153 . In both cases one arrives 
at essentially the same topological signature. 

Figure 2. The contribution $^(sj) of the isometrics g — a [part (a)] and g — /3 [part 
(b)] to the topological signature ^p'^{si) of the cubic manifold 71 with L = 1 in 
the observed universe with radius a ~ 0.71 . The horizontal axis gives the pair 
separation s while the vertical axis gives the normalized number of pairs. The 
MPSH approach was used to obtain both MPSH's <$^(si) > . These figures make 
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clear that the action of the isometrics a and (5 starts, respectively, at s = 1 and 
s ~ 0.71. These starting points for the actions of those isometrics give rise to the 
two discontinuities of the topological signatures ^^{^Si) of figures la and lb. 

Figure 3. The EPSH $e^p(sj) as given by eq. (p.22|) for an Euclidean universe Ba with 
radius a ^ 0.71 and fulfilled with an uniform distribution of cosmic objects. The 
horizontal axis gives the pair separation s while the vertical axis gives the normalized 
number of pairs. 

Figure 4. The contribution of the terms v^'"(sj) [part (a)] and V'^(sj) [part (b)] as given, 
respectively, by eqs. ( p.l7| ) [or (pI7|)] and ( |2.16| ), to the topological signature ^^{si) 
of the cubic manifold 7^ with L = 1 in the observed universe with radius a ~ 0.71 . 
The horizontal axis gives the pair separation s while the vertical axis gives the 
normalized number of pairs. The exact expression (p.22|) was used in part (a) and 
MPSH technique was used in both part (a) and part (b). 

Figure 5. The sum of the contributions due to the terms ^i^Si) and y9'^(sj), given in 
figure 4, to the topological signature v^'^(sj) [as given by eqs. (|2.15|) - ( p.l7| )] of 
the cubic manifold T\ with L = 1 in the observed universe with radius a ^ 0.71 . 
The horizontal axis gives the pair separation s while the vertical axis furnishes the 
normalized number of pairs. The comparison of this figure with figure 1 shows that 
the two procedures for obtaining the topological signatures give rise to the same 
topological signature. 

Figure 6. The topological signature ^^{si) = (n — 1) [< $(sj) > — $g^p(si)] for an elliptic 
manifold S^/Z^ with the edge of the tetrahedron L ~ 1.82. The observed universe 
was taken to be the whole unitary sphere S^, and the expression ( |2.24D for $g^p(sj) 
was used. The horizontal axis gives the pair separation s while the vertical axis 
gives the normalized number of pairs. 

Figure 7. The topological signature f^{si) = {<n> — 1) [<$(sj)> — ^Hpisi)] for the 
Seifert- Weber dodecahedral (hyperbolic) space with L = H = 1.992769 (edges L 
and height H), where the lengths are measured with the hyperbolic geometry (|1.1| ) 
with R = 1. The observed universe Ba with diameter 2a = 2.88 was used. The 
expression ( p.23| ) for $g^p(sj) was used. The horizontal axis gives the pair separation 
s while the vertical axis furnishes the normalized number of pairs. 
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